Quantum Hall Effect from the Topological Surface States of Strained Bulk HgTe 
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We report transport studies on a three dimensional, 70 nm thick HgTe layer, which is strained 
by epitaxial growth on a CdTe substrate. The strain induces a band gap in the otherwise semi- 
metallic HgTe, which thus becomes a three dimensional topological insulator. Contributions from 
residual bulk carriers to the transport properties of the gapped HgTe layer are negligible at mK 
temperatures. As a result, the sample exhibits a quantized Hall effect that results from the 2D 
single cone Dirac-like topological surface states. 



PACS numbers: 

The discovery of two (2D) and three dimensional (3D) 
topological insulators (TI)^ has generated stro ng ac- 
tivity in the condensed matter physics communit}^^. 
Current research on 3D TIs is mostly focused on Bi2Te3, 
Bi2Se3 and Sb2Te3 compound J^HI^ due to their simple 
Dirac-like surface states, which have been observed by 
angle resolved photoemission spectroscopy (ARPES) and 
scanning tunneling microscop}'^^. However, these com- 
pounds show strong defect doping and low carrier mo- 
bility, and the observation of surface charge transport 
is obscured by the bulk conductivity. Many of the pre- 
dicted novel properties of a 3D TI, such as the quan- 
tized magneto-electric effectP^^ and the surface Majo- 
rana fermionJI^, can only be observed when bulk carriers 
are negligible compared to the surface states. Experimen- 
tally reaching the intrinsic TI regime, where bulk carriers 
are absent, is now the central focus of the field. 

The two dimensional TI state was first predicted 
and observed in 2D HgTe quantum wells (QW)P^, 
and non-local transport measurements demonstrate 
edge state transport without any contributions from 
2D bulk carriers. 3D HgTe is a semi- metal which is 
charge-neutral when the Fermi energy is at the touching 
point between the light-hole and heavy-hole bands 
at the Brillouin zone center. A unique property of the 
band structure of HgTe is the energetic inversion of 
the Tq and Fg band ordering, which is the origin of 
the quantum spin Hall effect in 2D HgTe/CdTe QWs^. 
Due to the band inversion, 3D HgTe is also expected 
to have Dirac-like surface stateJ^l^, but since the 
material is semi- metallic, this state is always coupled to 
metallic bulk states. With applied strain, a gap opens 
up between the light-hole and heavy-hole bands , so 
that strained 3D HgTe is expected to be a 3D T&^. 
In this paper we demonstrate experimentally that a 
gap opens up in in-plane strained 3D HgTe bulk layers 
grown by molecular beam epitaxy (MBE), and we reach 
the much sought after intrinsic TI regime in a material 
with negligible bulk carriers. In this regime, the Hall 
effect of the 3D HgTe bulk layer is quantized, due to the 
contributions from the surface states only. Theoretical 



considerations are in agreement with the experimental 
results and confirm the transport through 2D surface 
states with Dirac type dispersion. 

HgTe bulk samples have been grown by MBE on CdTe 
subt rates, which have a lattice constant that is 0.3 % 
larger than that of bulk HgTe (0.646 nm). At this 
mismatch, the critical thickness for lattice relaxation is 
around 200 nm, implying that for thinner HgTe the epi- 
layer adopts the lateral lattice constant of the substrate, 
while in thicker layers the strain is relaxed by the forma- 
tion of dislocations. 

To provide evidence for the occurrence of the topolog- 
ical surface state, we first show an ARPES measurement 
on a 1 /im thick HgTe layer in [Fig. 1 (a)]. In this layer, 
the lattice strain is fully relaxed, and the surface has the 
lattice constant of unstrained bulk material. The figure 
clearly shows the presence of the predominantly linearly 
dispersing surface state band (SSB), coexisting with bulk 
bands (BBs) (more data on these assignments can be 
found in the supplementary material). According to the 
theoretical analysis in Ref. (19), the surface states orig- 
inate from the inversion between the Fg light-hole and 
Tq electron bands, while the bulk bands, which appear 
nearly at the same energy range with the surface states, 
correspond to the Fg heavy-hole band. 

Since the relaxed sample is a semi-metal and thus not 
a TI in the strict sense, we have grown a thinner sample, 
which is only 70 nm thick, thin enough for the epitaxial 
strain due to the lattice mismatch with the substrate 
to coherently strain the sample, thus opening up a bulk 
insulating gap. Fig. 1 (b) shows a high resolution X-ray 
diffraction map of the [115] reflex of this sample in units 
of reciprocal lattice space vectors. The bright spot in 
the center of this graph is the reflex from the substrate, 
while the thin vertical line stems from the HgTe epilayer. 
The absence of any relaxation of the reciprocal lattice 
vector Qx for the epilayer is direct evidence that this 
sample is fully strained. The deformation potentials of 
(Hg,Cd)Te have been reported in the literatur^^ using 
these values and the 0.3 % lattice mismatch we calculate 
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an energy gap of the order of ~ 22 meV for fully strained 
HgTe on CdTe, using an eight band k-p model^^. 

For transport experiments, we have subsequently pat- 
terned parts of both wafers into Hall-geometry devices 
with a mesa of 200 jam width and 600 jam length, using 
argon ion etching. The magneto transport of the samples 
has been investigated at a base temperature of 50 mK, 
in magnetic fields up to 16 T. In the 1 /im thick sample, 
the Hall data indicates that bulk conductance dominates 
the transport. Fig. 1 (c) shows that for this sample we 
observe a non-monotonic dependence of the Hall voltage, 
which is characteristic for the multi-carrier transport ex- 
pected from a semi-metal. Much more interesting behav- 
ior is observed for the 70 nm-thick sample, of which the 
longitudinal and Hall resistance are shown in the inset 
of Fig. 2. From the low field data, the electron mobil- 
ity can be extracted, and estimated as 34000 cm^/(V-s), 
which is significantly higher than that observed in Bi2Se3 
and Bi2Te3^^'^. At high magnetic fields, the longitudi- 
nal and Hall resistance exhibit distinct features which 
are characteristic for a 2D electron system (2DES): the 
Hall resistance Rxy shows plateaus at the same magnetic 
fields where the longitudinal resistance Rxx develops a 
minimum (inset of Fig. 2). Additionally, the Hall resis- 
tance Rxy shows the expected 2D quantized plateau val- 
ues, which become clearer in a conductivity plot (Fig. 2). 

Compared with the quantum Hall effect of an ordinary 
2DES, two unusual observations should be emphasized, 
which indicate that the observed quantum Hall plateaus 
indeed result from the Dirac-type dispersion of the topo- 
logical surface states. First, at low magnetic fields, a 
sequence of Hall plateaus develops with odd filling fac- 
tors, = 9, 7, and 5, before at higher field the sequence 
is continued with z^ = 4, 3, and 2. The occurrence of 
odd-number filling factors at low magnetic field indicates 
the presence of a zero mode (a Landau level at zero en- 
ergy) due to the linear dispersion of Dirac fermions, as 
has also been found in graphend^ and in HgTe/CdTe 
QWs with a critical thickness of 6.3 nnP^. A 70nm thick 
layer can safely be regarded as a 3D material; the con- 
finement energies of bulk carriers are sufficiently small 
that Hall quantization effects would be washed out by 
multi-subband averaging. We thus assume that the Hall 
plateaus result from the topological surface states, and 
use a model with two Dirac cones, one on the top (vac- 
uum) and one on the bottom (CdTe interface) to de- 
scribe the system (see additional online material for de- 
tails). Because massless Dirac fermions in a magnetic 
field always exhibit a zero mode, the Hall conductance of 
a single Dirac cone is given by axy = (n + ^)e^//i, where 
n = fit OT n = rib (for the top and bottom surfaces, re- 
spectively) is always an integei^l^. The fractional factor 
of ^ in the Hall conductance axy is a consequence of the 
quantized bulk topological term in the electromagnetic 
action^^, and is independent of the microscopic details. 
When top and the bottom surface have the same filling 
factor, i.e. n = = n^, the total Hall conductance is 
given by <jj^ = (2n + l)e^//i. Therefore, within the two 



Dirac cone model, the odd filling factor at low magnetic 
fields can be naturally explained assuming both surfaces 
have the same density. The appearance of an even filling 
factor at high magnetic field indicates that the degener- 
acy is removed. In an ordinary 2D electron gas, such a lift 
of degeneracy may occur due to Zeeman coupling. How- 
ever, if a system with top and bottom surface states has 
inversion symmetry, Zeeman coupling cannot split the 
energies of the two remote Dirac cones. As we explain in 
the supplementary material, this is because the inversion 
symmetry is preserved by the magnetic field. The Lan- 
dau levels from the top and bottom surface states will 
remain degenerate as long as the hybridization between 
the two surface states is negligible. For a thickness of 
70nm, hybridization between the top and bottom sur- 
faces can be neglected since the surface state width is 
around 2~3nm. Thus we conclude that the inversion 
symmetry breaking is necessary for the explanation of 
the Landau level splitting. In practice, the different elec- 
trostatic environments of both surfaces break the inver- 
sion symmetry and lead to different densities at both 
surfaces. This then leads to different Landau filling fac- 
tors for the top and bottom surfaces at high magnetic 
fields, and the visibility of even filling factors at high 
fields results from the increased energy splitting between 
adjacent Landau levels at lower filling factors in a Dirac 
system, as schematically indicated in the inset of Fig. 
3. A calculation of the density of states in a magnetic 
field from a two Dirac cones model with an inversion 
breaking term indeed agrees well with the minima of the 
Shubnikov-de Haas oscillation, as shown in Fig. 3 (see 
additional online material for details). From the fit, we 
obtain very reasonable carrier densities of 3.7 x lO^^cm"^ 
for the bottom (CdTe interface) and 4.8 x lO^^cm"^ for 
the top surfaces, respectively. We note that the carrier 
densities found here also imply we can neglect any ef- 
fects from bulk carriers - if the amount carriers found 
here would be evenly distributed over the 70 nm slab 
this would result in an overall 2D density of ~ lO^^cm"^, 
yielding a very different quantum Hall behavior. 

The second observation is that the minima in Rxx ^o 
not approach zero even at the highest magnetic fields al- 
though the Hall resistance is quantized and the plateaus 
have the expected resistance value. This indicates that 
besides the chiral edge modes from the quantum Hall 
effect, there are other modes contributing to the longitu- 
dinal, but not the Hall transport. A plausible candidate 
are the metallic states at the side surfaced. While the 
above two Dirac cone model only takes the top and bot- 
tom surfaces explicitly into account, the topological sur- 
face states also occur on the four side surfaces, which see a 
parallel, rather than perpendicular, magnetic field. Con- 
sequently, at the side surfaces the Dirac points are not 
gapped, but only shifted by the applied magnetic field. 
Thus, the surface states at the side surfaces will remain 
metallic in magnetic field and coexist with the chiral edge 
states, as shown schematically in Fig 1(d). This provides 
a backscattering mechanism when the transport on the 
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FIG. 1: Supporting characterization experiments, (a) 

ARPES measurements on a relaxed, 1 /xm thick HgTe sam- 
ple. The dispersion of the surface state (SSB) and bulk bands 
(BB) are indicated by the arrows; (b) Reciprocal space map 
in the region of the [115] reflex of the 70 nm thick HgTe sam- 
ple demonstrating that the epilayer is coherently strained; (c) 
Hall resistance as a function of magnetic field of the semi- 
metallic 1 /im thick sample ; (d) Schematic picture of the 
coexistence of the chiral edge states from the upper and lower 
surfaces and the non-chiral metallic modes at the side sur- 
faces. The magnetic field is perpendicular to the upper and 
lower surfaces, but parallel to the four side surfaces. 
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FIG. 2: Transport data on the strain-gapped 70 nm 
thick HgTe sample. The Hall conductivity of the 70 nm 
thick HgTe sample measured at 50 mK shows plateaus at the 
quantized values. The inset shows the Hall resistance Rxy, 
together with the longitudinal resistance Rxx- 



Appendix A: Details of the ARPES Experiment 



top and bottom surfaces is in the quantum Hall regim^^. 
A more systematic study of the influence of the residual 
bulk conductivity and side surface states is required to 
fully understand the quantitative behavior of Rxx and 
Rxy, which is beyond the scope of the present letter. 

In summary, we have experimentally reached the much 
sought after regime of the intrinsic 3D topological insu- 
lator with negligible bulk carriers, in epitaxially strained 
3D HgTe sample. Our observation of a quantized Hall 
conductance in a 3D sample conclusively demonstrates 
a key feature of 3D topological insulators. A simple 
model with two Dirac cones is proposed to understand 
the most salient features of the transport measurement 
qualitatively. The quality of our sample should be suf- 
ficient to observe the quantized topological magneto- 
electric effect and directly determine the 3D topologi- 
cal invariant experimentally. 

We acknowledge useful discussions with G. Astakhov, 
J.E. Moore and A.H. MacDonald. This work was sup- 
ported by the German Research Foundation DFG (SPP 
1285 Halbleiter Spintronik, DFG-JST joint research pro- 
gram, and grants AS327/2-2 (E.G.N.) and HA5893/1- 
1(E.M.H.)), the Alexander von Humboldt Foundation 
(C.X.L. and S.C.Z.), the EU ERG-AG program (L.W.M.) 
and the US Department of Energy, Office of Basic Energy 
Sciences, Division of Materials Sciences and Engineering, 
under contract DE-AC02-76SF00515 (S.G.Z). 



Angle resolved photoemission spectroscopy (ARPES) 
experiments on 1/im thick HgTe thin film (grown on 
GdTe substrate) were performed at Beamline 10.0.1 of 
the Advanced Light Source (ALS). During the exper- 
iments, the measurement chamber pressure was kept 
< 4 X 10~^^ Torr, and data were recorded by a Scienta 
R4000 analyzers at 15 K sample temperature. The pho- 
ton energy used was 65 eV and the total convolved en- 
ergy and angle resolutions were 20 meV and 0.2 degree, 
respectively. The fresh surface for ARPES measurement 
was obtained by gently sputtering the HgTe (001) surface 
with Ar ion-beam, and no observable surface degrada- 
tion was noticed during typical experimental period (12 
hrs). We also measured the photon energy dependent 
(65-70eV) band dispersion of HgTe, as shown in Fig [4] 
For all measurements at different energies, the surface 
state band can clearly be observed, superimposed on a 
blurry background of bulk bands. The shape and disper- 
sion of the band do not vary with photon energy, while 
its intensity may vary. This which is direct evidence that 
the feature does not have dispersion, as expected from 
a two dimensional (surface) state. We note that at the 
vicinity of the Gamma point, the surface state band has 
gained finite curvature, similar to that observed^ in 
Bi2Se3. The difference between the HgTe and the Bi2Se3 
cases is that for the latter, both the upper and lower 
Dirac cone can be seen, while for the former case, the 
Dirac point resides close to the Fermi energy and thus 
only the lower Dirac cone can seen in the ARPES mea- 
surements of Fig. |4j 
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FIG. 3: Comparison between the calculated density 
of states (DOS) and the measured Shubnikov- de 
Haas oscillations. The DOS are calculated for a two Dirac 
cones model for the surface states with one Dirac cone (at 
the CdTe interface) having a carrier density 3.7 x lO^^cm"^ 
(blue dashed-dotted line) and the other (the surface facing 
vacuum) 4.8 x lO^^cm"^ (green dashed line). The sum of the 
DOS for the two Dirac cones (red line) compares well with 
the measured longitudinal resistance Rxx (black line). The 
inset shows the Landau ladders of the Dirac fermions on both 
surfaces. 



Appendix B: k«p calculation of the band structure 
of strained HgTe 

To understand the transport properties of the 70nm 
sample, we first calculate the band dispersion of the 
structure using the k • p type approach. HgTe can be 
well described by the eight-band Kane model, where the 
electron Tq band and hole Ty, Tg bands are taken into ac- 
count explicitly. The detailed form of the Kane model, as 
well as the related band parameters, can be found in Ref 
[22l As is evident from the reciprocal lattice map of Fig. 
lb in the main article, the 70 nm thick HgTe epilayer is 
fully strained due to the lattice mismatch between CdTe 
substrate and HgTe. Although this lattice mismatch is 
small (e = {acdTe - ciHgTe) / dRgTe ~ 0.003 with acdTe 
and a HgTe lattice constants of CdTe and HgTe, respec- 
tively), the resulting strain is sufficient to open a bandgap 
between the heavy- and light-hole bands. The CdTe sub- 
strate is oriented along (001), and for this orientation the 
strain tensor components are given by Cxx = ^yy = 
Czz = — 2eCi2/Cii and e^^j = for = x^y^z. Here, 
Cii and Ci2 are elastic stiffness constants. Because the 
off-diagonal components of the strain tensor (the shear 
components) are zero, there are no internal electric fields 
generated in the HgTe layer due to piezoelectric effects. 
The strain tensor components for an arbitrary growth di- 
rection can be determined using the model of De Caro et 
aP^. The effects of the strain tensor can be incorporated 
in the Kane model through the Bir-Pikus Hamiltonian^^, 
which can be easily obtained from the Kane Hamilto- 
nian with the substitution kikj Cij. For the present 
case, one finds that lattice strain shifts the light-hole 



and heavy-hole band-edges at the F-point and leads to 
the opening of a gap between them which is given by 
Eg = —b{exx + ^yy — '^^zz)^ whcrc b is the uniaxial de- 
formation potential, which amounts to -1.5V for HgT^^. 
For C12/C11 = 0.68^^ we obtain Eg ^ 22 meV. Using 
the above model, we have calculated the band structure 
of a 0.3 % strained 70nm thick HgTe slab, yielding the 
dispersion shown in Fig. |5j At F point, the conduction 
band edge is around llmeV while the valence band edge 
around -llmeV, and the direct gap at F point amounts to 
22meV. Within the bulk gap, we find two special states 
( shown in blue in Fig.l) which touch the valence band 
states at the F point, but merge with the conduction 
band for large k. In the present system, the Dirac point 
of the surface state is buried deep within the heavy hole 
valence band, about 70meV below the valence band edge 
for the HgTe-CdTe interface and lOOmeV for the HgTe- 
Vacuum surface, as shown schematically by the red guid- 
ing lines in Fig. [sj Therefore, these special states (shown 
in blue in Fig.l) consist of the surface states originating 
from the Dirac type dispersion between Tq electron and 
Fg li ght-ho le bands hybridized with the Fg heavy hole 
bandPEH. 

In actual devices, positions of the Dirac points can 
vary from ones shown here. This is due to band bending 
from the Hartree potential, which has not been taken into 
account in the present calculations. Obviously, the exis- 
tence of the surface states is not sensitive to the Hartree 
potential due to their topological nature. 
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FIG. 4: a-e, Energy dependent band structure along M-F- 
M direction from 65-70eV excitation photons, with the Ef, 
surface state band (SSB) and the bulk band (BB) marked. 



Appendix C: Two Dirac cone model 

In this section, we will introduce an effective model 
with two Dirac cones to describe the surface states of a 
slab of a topological insulator and investigate the Lan- 
dau level structure evolving from such a system when a 
perpendicular magnetic field is applied (this defines the 
z-direction). On the surface of a strong topological in- 
sulator, the low-energy physics of the surface states can 
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be described by the Dirac Hamiltonian, with two com- 
ponents related to each other by time reversal, which we 
refer to as spin for simplicity. In our slab, we need to 
take two surfaces, which are perpendicular to the mag- 
netic field and related to each other by inversion, into 
account, and consequently, there are in total two Dirac 
cones and the four basis states can be written as |a,cr), 
where a = ± denotes the upper and lower surface and 
(J =t,i denotes spin. Within this four state basis, the 
effective Hamiltonian is written as 

= hvfikxCTy - kyCFx) (g) r^, (CI) 

where the Pauli matrix a describes spin and r the op- 
posite surfaces; the Hamiltonian is invariant under the 
inversion operation P = 1 r^^ and the time reversal op- 
eration T = —icFy ^ where K indicates complex con- 
jugation. We note that for zero gap HgTe/CdTe quantum 
wellJ^, the effective Hamiltonian can also be described 
by two Dirac cones. However in that case, the two Dirac 
cones are related to each other by time reversal, while 
here the two Dirac cones are related by inversion. Be- 
sides the effective Hamiltonian H2D^ we also take into 
account two more terms that may be active in the sam- 
ple: an inversion breaking term^^ 

Hib=(^^Q =Ail$5r„ (C2) 

and a hybridization term between the two surface states, 
which should become important in very thin slabs: 

Hh=(^^^ ^'') =^hl®T,.. (C3) 

To study the effect of a magnetic field along the z- 
direction, we first need to make a Peierls substitution, i.e. 
we replace k in the Hamiltonian H2D with tt = k + | A 
, where A = (0,Po^,0) for a magnetic field B = Bqz. 
Additionally, we need to consider a Zeeman type term, 
given by 

where is an effective g factor that takes orbital effects 
of nearby bands into account. 

First let us have a look at the Landau levels of the 
Hamiltonian H2d^ which can be solved as 



Eo,t{n) =t^2nehv)BQ n = l,2, ••• (C5) 

and Ecyin = 0) = for zero modes, where t = ± denotes 
the electron and hole levels and a = ± denotes the up- 
per and lower surfaces. Obviously, each Landau level is 
doubly degenerate since the the upper and lower surfaces 
are identical and related by inversion. The zero modes 



will be half filled at charge neutrality, and, thus due to 
the double degeneracy, the Hall plateaus should appear 
only at odd number filling factors n = l,3,5,---, which, 
as described in the main article, is what we observe in 
low magnetic fields. 




k, nm" 

FIG. 5: Band structure of a 70nm thick 0.3 % strained HgTe 
slab. The Dirac- like states in the gap are plotted in blue. 
The dashed red line schematically shows the dispersion of the 
Dirac surface states at the two opposite surfaces before their 
hybridization with the Fg heavy hole bands. 

In high magnetic fields. Hall plateaus at even num- 
ber filling factors are also observed, which indicates a 
splitting of the Landau level degeneracy. In zero gap 
HgTe/CdTe quantum wells (QWs) case^^, we have found 
that the even number filling factors observed in that sys- 
tem are caused by Zeeman splitting. It is thus natural to 
first investigate the same mechanism for the present ex- 
periment. This implies solving the Landau level problem 
for the Hamiltonian H2D + + Hz- One easily finds 
that in this case the Landau levels are given by 

Eat{n) = t^2nehv}Bo + (l/iB^^ol + a\Ah\)^ (C6) 

for n = 1,2,3, ••• and Ea{0) = -g* /ibBq a\Ah\ for 
n = 0. From these expressions, we conclude that the 
double degeneracy remains when only the Zeeman term 
is important but hybridization between the two surface 
states is negligible ( A /j, = 0). Therefore considerable hy- 
bridization between the two surface states (i.e. the slab 
has to be very thin) is required in order to make this 
mechanism work. In HgTe, the decay length of the sur- 
face states is estimated to be of order 2 ~ 3nm^^, which 
is much smaller than the sample thickness of 70 nm. It 
is thus unlikely that a Zeeman mechanism is responsible 
for the degeneracy lifting in the present experiment. 

The other possible mechanism for Landau level split- 
ting is the breaking of the inversion symmetry - which in 
practice occurs by the presence of the substrate. To see 
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FIG. 6: Landau level index for the data of fig. 2 in the main 
article plotted as a function of inverse magnetic field. 



this, we solve the Landau level problem of the Hamilto- 
nian H2D -\- and obtain the eigen-energies as 



Eat = ty^/2nehv'jBo + aA. 



(C7) 



for n = 1,2,3, ••• and Ea{0) = aAi for n = 0. From 
this expression, it is obvious that the degeneracy of the 
Landau levels is removed for any nonzero A^. This makes 
inversion-symmetry breaking mechanism the most likely 



explanation for our observation of even filling factor at 
high magnetic fields. 



For a comparison with the experimentally obtained 
Shubnikov-de Haas oscillations, we calculate the den- 



sity of states for the energy dispersion (C7) for two 



Dirac cones with an inversion-symmetry breaking term. 
Taking hvf = 280meV • nm and = A.2meV ^ we 
find that the positions of the experimental Shubnikov- 
de Haas minima can be fitted well with the densities of 
m = 3.7 X lO-^-^cm"^ and 712 = 4.8 x 10^^cm~'^ for the 
two Dirac cones, which are the traces shown in Fig. 3 
of the main article. Further evidence that this is the 
correct model for our observations comes from the odd 
filling factor {u = 9, 7, 5) quantum Hall plateaus we ob- 
serve at low magnetic fields. A clear discrimination be- 
tween the Landau level structure of a Dirac system and 
that of a conventional 2D electron gas can be obtained 
fro m plo tting of the Landau level index as a function of 
For a Dirac system with multiple Dirac cones, 
the filling factor is related to the Landau level index N 
hy u = m{N + 1/2), where m is the number of the Dirac 
cones and equals to 2 here. Taking the magnetic field val- 
ues corresponding to the Hall plateaus u = 9, 7, 5 from 
the axy curve in Fig. 2 in the main article, and plotting 
N as a function of 1/B yields Fig. [gJ The intercept of this 
plot for infinite magnetic field gives -1/2, which provides 
additional evidence that the main physics of our system 
can be well described by the two Dirac cones model. 
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